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ABSTRACT 

Aims. The present paper addresses some of the problems in the buildup of element stratification in stellar magnetic atmospheres due to 
microscopic diffusion, in particular the redistribution of momentum among the various ionisation stages of a given element and the calculation 
of diffusion velocities in the presence of inclined magnetic fields. 

Methods. We have considerably modified and extended our CARAT code to provide radiative accelerations, not only from bound-bound 
but also from bound-free transitions. In addition, our code now computes ionisation and recombination rates, both radiative and collisional. 
These rates are used in calculating the redistribution of momentum among the various ionisation stages of the chemical elements. A careful 
comparison shows that the two different theoretical approaches to redistribution that are presently available lead to widely discrepant results 
for some chemical elements, especially in the magnetic case. In the absence of a fully satisfactory theory of redistribution, we propose to use 
the geometrical mean of the radiative accelerations from both methods. 

Results. Diffusion velocities have been calculated for 28 chemical elements in a r eff = 12000 AT, logg = 4.00 stellar magnetic atmosphere 
with solar abundances. Velocities and resulting element fluxes in magnetic fields are discussed; rates of abundance changes are analysed for 
systematic trends with field strength and field direction. Special consideration is given to the Si case and our results are confronted in detail 
with well-known results derived more than two decades ago. 
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1. Introduction 

Having been proposed more than 3 decades ago by Michaud 
( 1970 1, the diffusion theory still constitutes the only viable sce- 
nario for the buildup of abundance inhomogeneities in chemi- 
cally peculiar (CP) early type main-sequence stars. Early cal- 
culations by Michaud et al. J1976i have shown that diffusion 
in the envelope of main-sequence stars can lead to overabun- 
dances of a number of chemical elements by up to 7 dex, which 
would be sufficient to explain even the largest overabundances 
found in CP stars. 

Magnetic fields can be expected to modify diffusion, and 
indeed Vauclair et al. 0979) outlined a particular mechanism 
involving upward diffusion of Sii that produces overabun- 
dances in magnetic stars at places where the magnetic field is 
horizontal. Whether vertical or horizontal magnetic fields could 
lead to the development of spots or rings with enhanced ele- 
mental abundances depends on the element in question; silicon 
appears to concentrate in horizontal field regions (Alecian & 



Vaucrah Tl981l l. whereas gallium is thought to accumulate near 
vertical field regions (Alecian & Artru 1987 ). 

The detailed modelling of 53 Cam by Babel & Michaud 
J1991bl has shown that a simple diffusion model is not al- 
ways sufficient and that other processes should be considered, 
like inhomogeneous mass-loss. Whenever one invokes diffu- 
sion of chemical elements as the mechanism responsible for 
abundance anomalies in a star, this implies that a stratifica- 
tion process has been at work somewhere in the star or is still 
at work, depending on the stellar age and its fundamental pa- 
rameters. Observational evidence for such vertical abundance 
stratification of chemical elements has accumulated over recent 
years (see Ryabchikova et al. 2002). Zeeman-Doppler-Imaging 



in all 4 Stokes parameters (Kochukhov et al. I2004t is begin- 
ning to provide simultaneous magnetic and abundance maps 
for several elements, making it possible to assess whether there 
is a correlation between the respective topologies of the surface 
magnetic field and of the chemical abundances, as predicted by 
diffusion theory. 
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As far as stellar interiors are concerned, the last decade has 
seen radiative acceleration calculations based on large opacity 
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and atomic databases, resulting in increasingly realistic mod- 
elling of the diffusion of chemical elements. In the outer layers 
of the stars, the treatment of radiative transfer proves far more 
difficult and time-consuming; likewise, large-scale calculations 
for a fair number of chemical elements, involving reasonably 
realistic spectra, the correct treatment of blends, and accurate 
formal solvers require very fast, preferably parallel computers, 
especially if Zeeman splitting and magnetic polarisation is to 
be taken into account (Alecian & Stift l2002l 120041 ) . Previously, 
Zeeman splitting was estimated to be of merely secondary im- 
portance in the case of Si (Alecian & Vauclair 1981 ), so only 
the effects of the magnetic field on the movements of the ions 
across magnetic lines were being considered (see Vauclair et al. 
119791 . Later Babel & Michaud ( 1991ai found moderate ampli- 
fications of the radiative accelerations (due to Zeeman splitting) 
over the zero field value. Alecian & Stift (2004i investigated 
the magnetic amplifications of a large number of chemical el- 
ements as a function of field strength and of field direction in 
detail. Radiative accelerations were shown to be quite sensitive 
to blending, magneto-optical effects, and the Zeeman patterns 
of the lines involved. 

Since our previous study of radiative accelerations in mag- 
netic fields (Alecian & Stift 2004) was mainly concerned with 
the amplification effect due to Zeeman splitting, we only con- 
sidered the radiative acceleration due to bound-bound transi- 
tions. In the present paper, we want to address the problem of 
element stratification in stellar magnetic atmospheres due to 
microscopic diffusion. Therefore, we have to compute as accu- 
rately as possible all the physical ingredients that determine the 
fluxes of the various particles and that eventually lead to abun- 
dance inhomogeneities. In the following section we present 
detailed computations of radiative accelerations, including the 
contributions of bound-free transitions and the redistribution 
of momentum among the various states of ionisation of a given 
element. In the subsequent section, we consider the diffusion 
velocity in the presence of inclined magnetic fields and discuss 
the basic mechanism for the build-up of abundance stratifica- 
tion. Finally, we outline the trends of abundance stratification 
in the magnetic atmosphere of a main-sequence star with effec- 
tive temperature T e g = 12000K. 



2. Radiative accelerations 

In Alecian & Stift ( 2004 1, we only considered the average ac- 
celeration of each element due to bound-bound transitions of 
ions. Actually, to compute the diffusion velocity, one also has 
to consider the acceleration due to bound-free interactions and 
correct the average acceleration for the redistribution of mo- 
mentum among the various ionisation states. 

2.1. Accelerations due to b-b and b-f transitions 

An ion A +l can be radiatively accelerated both through bound- 
bound transitions - discussed in detail in Alecian & Stift ( 2002 
2004 1 - and through bound-free transitions (photoionisations) 



of the ion A +(!_1) (see Alecian lT^t 

Xoo 
CTi,km <E>v dv (1) 

(1 -yi-u)vi-i,k® v dv (2) 

..-it. 

where a,- = l/^ m^ c). The mass of the atom is m^, the re- 
spective total and k-level number densities of ion A +l are de- 
noted by Af, and N^, and c is the velocity of light. Q> v is the 
radiative flux (in ergcirT 2 s Hz -1 ), cr, t m the cross-section of 
transition k — > m, and cr, the photoionisation cross-section 
with threshold frequency v,^. The correcting factory,^ has been 
introduced by Michaud J1970i and gives the fraction of mo- 
mentum imparted to the electron. Usually it is approximated by 
1.6 (1 ~Xi,k/hv) with Xi,k the ionisation energy (Sommerfeld & 
Schur. [T930l . Although recent analytical studies by Massacrier 
( 1996) have improved the knowledge of for H-like ions, it 
is poorly known for the rest, and we kept the usual approxima- 
tion. The fraction of momentum transferred to the ion is thus 1 
at threshold, goes to zero at hv = %Xi,k, and becomes negative 
at higher frequencies. There is no special treatment for autoion- 
isation transitions in our approach; some of them are included 
in the photoionisation cross-sections provided by TOPbase (see 
Sect.EJ. 

The CARAT code (for a description see Alecian & Stift 
2002 2004 1 carries out full opacity sampling of Zeeman split 
spectral lines and determines the radiative flux by solving the 
polarised radiative transfer equation. The model atmospheres 
used in our study have been calculated with the Atlas9 model 
atmosphere code (Kurucz 1993), so the continuous opacities 
are based on the Atlas9 opacity routines. The photoionisation 
cross-sections cr, ^ , however, are taken mainly from TOPbase 
as discussed below in Sect. [3] This introduces a slight incon- 
sistency: in contrast to the bound-bound case (Eq. Q) where 
cross-sections, fluxes, and resulting accelerations are based on 
the same set of atomic data - in our study we used data from 
VALD (Piskunov et al. 1995) - photoionisation cross-sections 
<x, a and flux Q> Y respectively in Eq. involve different atomic 
data sets. Still, we consider that such an inconsistency has 
no significant consequences on the diffusion velocities in the 
present work, since Atlas9 models have proven to yield fairly 
realistic stellar fluxes. Random errors may, however, stem from 
uncertainties in the computed wavelengths of photoionisation 
thresholds, which can be much larger than the typical Doppler 
widths of spectral lines. Pending the completion of a new stellar 
atmosphere code currently being undertaken by Stift and col- 
laborators, we are convinced that at present there is no viable 
alternative to our approach. 

2.2. Redistribution of momentum among ions 

In Alecian & Stift ( 2004 ), we were mainly interested in Zeeman 
amplification and we approximated the acceleration of a given 
element by the simple average value of the accelerations of its 
various ions, i.e. by the sum of the individual accelerations 
g md of the ions, weighted by their respective relative popula- 
tions. This was justified in that particular context, since we only 
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wanted to estimate the effect of magnetic fields on the acceler- 
ations. However, this kind of simple average is based on an 
implicit hypothesis, viz. it supposes that from the point of view 
of the transport process, the ions of a given element behave like 
particles that are completely independent of each other. In other 
words, it is assumed that they interact only by means of elas- 
tic collisions. Ionisation/recombination processes, however, are 
not elastic interactions, in particular because particles are not of 
the same type before and after the interaction. What is called 
population A +l actually designates (in the statistical sense) the 
average fraction of atoms A that are, at a given time, in the 
ith state of ionisation. From the microscopic point of view, in a 
stellar layer where ion A +l is dominant, the atoms A spend most 
of their life in state i, but also to a lesser degree in the adjacent 
ionisation states i - 1 and i + 1 . If momentum is gained by an 
ion in state i through absorption of a photon, and if the time re- 
quired to lose this momentum (i.e. transfer it to the medium by 
means of collisions) is larger than the mean time for ionisation 
or recombination, a non-negligible part of this momentum will 
be shared with states i — 1 or i + 1 . This is approximately equiv- 
alent to the assumption that the corresponding atomic radiative 
force will be exerted mainly on state i but also to a lesser de- 
gree on the ionisation states i — 1 or i + 1, which result from 
recombination and ionisation respectively. 

This picture applies to the momentum gained through 
bound-bound transitions. In a bound-free transition, one as- 
sumes that the ionisation time is zero, so most of the momen- 
tum is acquired by ion A' +1 but for the momentum carried away 
by the ejected electron. That immediate momentum redistribu- 
tion for bound-free photoabsorption is implicit in Eq.(0- 

The redistribution effect was first introduced by Montmerle 
& Michaud (119761 hereafter referred to as MM). It must be 
taken into account in the diffusion process, especially when the 
neutral state of an element is involved, because the diffusion 
coefficient for neutrals is several orders of magnitude larger 
than for ions (in layers where hydrogen is strongly ionised). 
Actually, redistribution of momentum raises complex physical 
questions that have not yet been completely solved. One of the 
main difficulties lies in the fact that the redistribution of mo- 
mentum must depend on the energy levels, since the probabil- 
ity of ionisation is higher from higher levels than from lower 
ones. In the original work of MM, this was not considered. 
Gonzalez et al. ( 1995 hereafter referred to as GLAM) studied 
this problem by means of a detailed study of carbon interaction 
rates (in the optically thick case) as a function of the energy 
levels. These authors also proposed an approximate method of 
estimating the redistribution of momentum. They found that 
the momentum of Cn is mainly redistributed to Cm when the 
absorption of photons leads to final energy levels with princi- 
pal quantum numbers greater than or equal to 3, which cor- 
responds to levels with an energy that is higher than 2/3 of 
the ionisation potential, and proposed to generalise these find- 
ings to all elements. Hereafter, we shall call this method the 
light-GLAM method. Proffitt et al. J1999I applied the detailed 
GLAM method to the study of Hg throughout the atmosphere. 
This is certainly the most accurate calculation of the radiative 
acceleration ever done for stellar atmospheres since the de- 
tailed interaction rates have been calculated for each energy 



level of Hg. In the present work and according to the present 
status of our CARAT code, we cannot apply such a sophisti- 
cated method. We therefore computed the radiative accelera- 
tions using both the MM and the light-GLAM method (see the 
discussion in Sect. 15. U . 

Details of the formulae, and equations used for estimating 
the respective interaction rates (or time scales) - collision, ion- 
isation, and recombination - and for computing the radiative 
accelerations with both the MM and the light-GLAM approach 
are presented in Appendixes|A]|B] andlCl 

3. CARAT : added input physics 

The CARAT code developed by Stift (Alecian & Stift |2"UU21 
2004) and based on the COSSAM code of Stift dSUUDl l. also 
discussed in Wade et al. (2001 1, has in the past only been 
capable of calculating accelerations due to bound-bound 
transitions. Since bound-free processes had not been included, 
actual diffusion velocities could not be determined. For the 
present study, we have therefore modified and extended 
CARAT to include routines that calculate photoionisation and 
collisional ionisation rates (see Appendix [Cji using several 
different sources. Tables of photoionisation cross sections 
as a function of energy come from TOPbase (Cunto & 
Mendoza 119921 Cunto et al. 1199 3): sufficient resolution of 
the resonances can lead to large sets of data points for many 
levels. Energy level data that can be used in conjunction with 
hydrogenic formulae (see also Alecian [1994) are extracted 
from TOPbase, from the NIST Atomic Spectra Database 
( http://physics.nist.gov/PhysRefData/ASD/index.htmli and 
from Kurucz ( 1993 ). 

As pointed out above, TOPbase data for ions A +l with 
charge i (http://vizier.u-strasbg.fr/TOPbasei consist of either 
detailed photoionisation cross section tables with 

- kS LP the quantum numbers of the spectroscopic series 
containing the level k, coded as (2S + 1)* 100 + L* 10 + P 
{P = 1 corresponds to odd parity), 

- kLV the energy position of level k within its spectroscopic 
series, 

- Er the level energy in Rydbergs (relative to the ionisation 
limit), 

- NP the number of subsequently listed energy - cross- 
section pairs, 

or of level listings which, in addition to kS LP, kLV and Er, 
include 

- gij,, the statistical weight of level k, 

- the effective quantum number. 

In contrast, level data from NIST or Kurucz give 

- J the total quantum number of level k, 

- E g the level energy (in cm -1 ), measured from the ground 
state. 

CARAT does not require any of the further atomic level data 
that are provided by TOPbase, NIST or Kurucz, such as nuclear 
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charge, number of electrons, electron configuration, spectro- 
scopic term, quantum defect, Lande factors, and radiative life- 
times of the levels. Statistical weights of the levels, when not 
listed, are either calculated from kSLP or from /. 

When only energy level data are available, we resort to a 
hydrogenic approximation of the photoionisation cross section 
(in cirT 2 ) from level k 



<r ik =7.910" 



»q tXi,k \ 3 
i + l) 2 I hv I 



(3) 



where i is the charge (the degree of ionisation) of the atom with 
i = for neutrals, Xi,k the ionisation potential measured from 
level k, and hv the energy of the photon. 

CARAT offers several options for the inclusion of pho- 
toionisation cross section data from the various sources. 
Presently we always use the complete set of TOPbase level 
data. Whenever detailed cross section data are available, they 
replace the level data - for the purpose of identifying the lev- 
els to be replaced we need the quantities kSLP and kLV. The 
elements covered by TOPbase are Li, Be, B, C, N, O, Ne, Na, 
Mg, Al, Si, S, Ar, Ca, and Fe. The NIST data are taken solely 
for those remaining elements where data are available for at 
least the first 3 ionisation stages, viz. F, P, CI, K, Sc, Ti, V, Cr, 
Mn, Co, Ni, Cu, Zn, and Ga. Fe constitutes a special case in 
that TOPbase does not list the first 2 ionisation stages; they are 
taken from NIST. We have omitted the rare earths. 

4. Diffusion velocity and element stratification 

The diffusion velocity Vd, of an ion in a binary mixture (pro- 
tons, electrons, and trace ions), without magnetic field (Aller & 
Chapman 1960 Michaud 1970 1, can be approximated by: 



V D , * Di 



where 



Sine; 
dr 



Am„ , 



rad 



r) 



(Z; + l)m p 
2kf 



(4) 



- D, is the diffusion coefficient of the ion i, 

- A m v the mass (cgs) of the atom, 

- r,T,k are the geometrical depth, the temperature, and the 
Boltzmann constant, respectively, 

- c, is the ion concentration (the ratio of the ion partial pres- 
sure over the total gas pressure), 

- Z, the charge (Z, = for the neutral state) 

- gi the redistributed radiative acceleration of the ion, 

- g the gravity. 

The third term in the brackets, introduced by Aller & Chapman 
( 1960 1, constitutes a minor correction due to the effect of the 
electric field. In this work, we have neglected thermal diffu- 
sion; the first term in brackets - ordinary diffusion due to a 
gradient in the concentration of ions - is small (due to the ion- 
isation gradient only) since in the present study we are consid- 
ering homogeneous abundances. A more sophisticated method 
of computing microscopic diffusion for a multi-component gas 
has been used in stellar evolution codes (Richer et al.[T998). 
The approximate Eq. © assumes that ions are trace particles 
and is accurate enough for our purposes. The diffusion coeffi- 



and Michaud & Proffitt ( 1993); here we use the same approx- 
imations as Landstreet et al. (1998) and Budaj & Dworetsky 



( 2002 ), corresponding to the adoption of the coefficient for bi- 
nary ion-proton collisions, corrected for collisions with neutral 
hydrogen. The latter correction could prove significant in at- 
mospheres with effective temperatures lower than 10000K. 

In Eq. ©, we have also neglected ambipolar diffusion 
(Babel & Michaud. ll991a> . This effect might be important in 
places where 7Y(Hi)/7Y(Hn) * 1 ; however, this does not occur in 
the model we will consider. 

4.1. Diffusion velocity in magnetic fields 

The approximate theory of diffusion in magnetic fields can be 
found in Chapman & Cowling (1970). The diffusion velocity 
@ of charged particles, orthogonal to magnetic lines, is re- 
duced by the factor: 



fs, 



low,i 



(1 



(5) 



where f,- is the collision time (the time taken by a particle i with 
mass m, and charge Z,e, to deviate by | through collisions from 
its initial motion), w,-/27r the cyclotron frequency in a field with 
strength H: 



ZteH 



(6) 



The collision times f, can be estimated easily from the diffusion 
coefficients (see AppendixlAV 

The average vertical diffusion velocity (in the presence of 
horizontal magnetic fields) of an element can be approximated 
by the following expression (the sums are over the ionisation 
states /): 



Vn 



2 fslowj VDi 



(7) 



cients have been discussed in detail by Paquette et al. ( 1986 



where the Ni are the number densities of the respective ions. 

This formula is the one used by Vauclair et al. !W79\ to 
study the Si accumulation in magnetic Ap stars. It is worth- 
while to recall here how these authors have analysed the be- 
haviour of the diffusion velocity as a function of magnetic field 
direction, since it is very instructive. Silicon is found to be over- 
abundant in magnetic Ap stars but more or less normal in non- 
magnetic peculiar stars. Vauclair et al. ( 1979 ) have shown that 
the radiative acceleration on ionised Si is not very strong (due 
to saturation of lines) and cannot therefore lead to high over- 
abundances, even though the acceleration on the neutral state 
is strong. This explains the quasi-normal abundance of Si in 
the non-magnetic stars. But in magnetic Ap stars, the down- 
ward diffusion flux due to the settling of Si ions (mainly Sin) 
is slowed down by the magnetic field (small f s i ow j in Eq. (0), 
while for the neutral state, f S hw,o = 1 ensures a positive veloc- 
ity Vd and leads to Si overabundances. The authors concluded 
that the diffusion model is compatible with the Si abundances 
observed in peculiar stars. Their study has been extended to 
the oblique rotator model by Michaud et al. (1981), assum- 
ing a dipole plus quadrupole magnetic structure. These authors 
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showed how elements can accumulate, depending on the orien- 
tation of the magnetic field (at the magnetic poles and/or equa- 
tor). 

So far, our considerations have essentially been restricted 
to the average diffusion velocity across horizontal magnetic 
lines. However, at the surface of magnetic stars, the direction 
of the magnetic field lines varies from to | between magnetic 
pole(s) and equator(s). This angular dependence can explain in- 
homogeneous abundance distributions over the stellar surface: 
the sensitivity of the diffusion velocity to the inclination of the 
magnetic field lines determines the size of the abundance struc- 
tures (see Michaud et al. 119811 and Alecian & Vauclair 1981 1. 
In magnetic fields that are neither strictly horizontal nor strictly 
vertical, the diffusion velocity vector is now inclined with re- 
spect to the vertical (z-axis). The inclination is given by the 
expressions derived by Alecian & Vauclair ( 1981 1 (please note 
that their 9 angle was defined with respect to the horizontal and 
not to the vertical as in the present work). The vertical (z-axis) 
component is: 



V H: ,i ~ V Di [fslowj + (1 - f llow ,i) cos 2 e] . 
and the horizontal one (x-axis) is: 

V Hx ,i ~ ~ (1 - Lhw,d V Di sin 26*. 



(8) 



(9) 



The diffusion velocity in a magnetic field is given by the 
vector Vhj = (VhxJ, Vhzj)- F° r vertical magnetic lines (8 = 0), 
one verifies that V//, = (0, Vd : ) according to Eqs. © and 
In this case, the only effect of the magnetic field comes 
from the amplification of the vertical component of radiative 
acceleration due to Zeeman splitting (Alecian & Stift 2004 1. 
For horizontal magnetic field lines (6 — |), one finds Vm = 
(0, f slowly d,) as in Eq. @. For inclined magnetic lines, the hor- 
izontal component given by expression ij^Ji represents only a 
horizontal diffusion of ions (keep in mind that the horizontal 
component of the radiative acceleration is generally negligible, 
see Alecian & Stift 2004 1. However, this horizontal component 
of the diffusion velocity, which could be of the same order of 
magnitude as the vertical one, has a very small effect on the 
horizontal distribution of elements. This results from a combi- 
nation of time scales and of spatial scales. Indeed, in the case of 
silicon, for instance, Megessier ( 1984 1 has estimated that Si mi- 
gration from the pole to the equator would take some 1 7 years ! 
Only if the life time of the magnetic structures is of this order of 
magnitude, can an observable effect of horizontal diffusion be 
produced. Because of the much smaller vertical height scales, 
vertical diffusion leads to observable effects much earlier (by 4 
to 5 orders of magnitude). One can thus conclude that it appears 
rather unlikely that horizontal diffusion will have a significant 
effect. This means that the abundance inhomogeneities on mag- 
netic Ap stars (patches or rings) are mainly due to the angular 
dependence of the vertical component of the velocity vector as 
given by Eq. combined with some local mass-loss velocity. 
Therefore we shall stick hereafter to the one-dimensional scalar 
form of Eq. Q and just use Vh,zI instead of /,/ow,; Vm- All the 
computations assume zero mass-loss velocity. 



4.2. Abundance stratification 

One of the main goals of Ap star modelling is the simula- 
tion of the development in time of abundance stratifications. 
Unfortunately, this is still beyond our possibilities, the strat- 
ification process being non-linear and requiring both further 
theoretical (and software) developments and huge computing 
power. Still, a first guess of the kind of stratifications to be ex- 
pected can be obtained by considering the flux of the various 
diffusing elements at the beginning of the process (f = 0), when 
elements still have their normal abundances and with no com- 
peting processes considered. We discuss this point in more de- 
tail in Sect. |5] Such a study was first carried out by Alecian 
& Vauclair (1981 1 who neglected the Zeeman effect. In the 
present work, by neglecting the horizontal flux as discussed in 
the previous subsection, we compute the vertical diffusion flux 
Fd, at depth z (increasing towards the centre of the star) by 

*D(z) = XiVfVW (10) 

from which we determine, for each element considered, the lo- 
cal rate of abundance change per unit time, at time t — 0. This 
rate is ao = d t N/N, where N is the number density of element 
A, i.e. Yj Ni- Let N s i a b be the total number of particles A in a 
cylinder of height Az and section unity. One has: 

3tN_ ^ A,N slab = F D (z + Az)-F D (z) 
N 



NAz 



N s l a b 

which leads to 
a =N- 1 d z F D (z). 



5. Results and discussion 



(11) 



(12) 



We have computed the radiative accelerations and diffusion ve- 
locities of 28 elements (Be to Ga) that are generally observed 
in Ap stars, and we have enough atomic data at our disposal for 
most of these to obtain reasonably accurate results. The main- 
sequence stellar model we adopted comes from Kurucz ( 1993 1 
with T e ff = 12000 K, \ogg = 4.00, and solar abundances. 

5.1. Computation of the radiative accelerations 

Hereafter, we consider the total radiative acceleration on each 
element, computed by both the MM and the light-GLAM meth- 
ods, with 



„rad 



rod 



(13) 



XNiD HzJ 
where 

D Hz .i = Di [f slowJ + (1 - fsio W ,i) cos 2 e] . (14) 

When the total radiative acceleration (vertical component) 
equals gravity and if concentration gradients are neglected, the 
diffusion velocity of the element should become zero (accord- 
ing to Eqs. Q and dlOt ). This leads quite naturally to the defini- 
tion of an effective total acceleration g^f (taking into account 
the correction due to the electric field): 



eff 
tot 



Z Ni D Hz j[gf-{l- z -£)g\ 
ZNiD HzJ 



(15) 
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-4 -3-2-10 1 2 

lo g T 5000 

Fig. 1. Comparison of the Si accelerations obtained by the 
light-GLAM method and the MM method (Eq. JT5})- The loga- 
rithm of gcij^M /&mm * s shown vs. the logarithm of optical depth 
for a model with T eW = 12000 K, log g = 4.00. The dashed line 
pertains to Gauss, the solid line to 10000 Gauss. 

Note that the g r " d are the respective redistributed accelerations. 
The total diffusion flux dlOi can be expressed directly as a func- 
tion ofg( { y : 

F D (z)*N(D)^g e t f/, (16) 
where we have defined the average diffusion coefficient as: 
(D) = J]N i D Hz j/N. (17) 

According to the study by GLAM, one expects to have large 
discrepancies between the MM and the light-GLAM methods, 
and we confirm this trend. In Fig.^we present the logarithm 
of the ratio of the accelerations light-GLAM vs. MM for silicon 
(to be discussed in some detail in Sect. 15. 3> . The discrepancy is 
especially strong in the magnetic case and for layers with opti- 
cal depth smaller than unity. We carried out such a comparison 
for all the 28 elements investigated, and found that the discrep- 
ancy is much smaller for some of them, but larger for some oth- 
ers. The elements for which the discrepancies are smaller than 
0.3 dex in the non-magnetic case are: Be, B, Ca, Sc, Ti, V, Cr, 
Mn, Fe, Co, Ni, Cu. The discrepancies are generally larger in 
the magnetic case. We do not notice obvious systematic trends 
except that MM is larger than light-GLAM for noble gases (Ne, 
Ar) and that MM is smaller than light-GLAM for the alkali el- 
ements (Na, K), i.e. elements with a noble gas configuration in 
their first ionisation stage. This reveals the role of the distri- 
bution of energy levels, since levels are not considered in the 
same way in light-GLAM and in MM. Silicon appears to be 
a rather typical case concerning the discrepancy between the 
MM and the light-GLAM methods; we discuss it in more detail. 
Figure^ shows that in the non-magnetic case, the MM accel- 
eration of Si is smaller by 0.66 dex, which is the same order 
of magnitude as the discrepancy found by GLAM for carbon - 
their discrepancy, however, goes in the opposite direction, their 
MM acceleration being stronger. Notice that we cannot com- 
pare our results concerning carbon with those of GLAM since 



their computations apply to a star with T e ff = 8000 K and to 
deep layers (outside the validity domain of CARAT). To better 
understand the reasons for such large differences, we examined 
the contributions of the various terms entering both methods 
in detail, especially for the neutral state since it plays a ma- 
jor role in momentum redistribution as already mentioned in 
Sect. 12.21 The light-GLAM method assumes that ionisation is 
dominated by transitions to the continuum from high energy 
levels, because after absorption of a photon leading to a suffi- 
ciently high-lying energy level, collisional excitations prevent 
de-excitation. For the MM method, we have considered the to- 
tal direct ionisation rate from any level (see Ea. JC.5» . without 
taking the added path of excitation collisions before ionisation 
into account. The MM method thus gives too much weight to 
low energy levels. If one considers that the MM method under- 
estimates the redistribution of momentum from the neutral to 
the first ionisation state, one is led to conclude that it gives an 
upper limit of the radiative acceleration. On the other hand, the 
light-GLAM method is certainly less justified in the upper at- 
mosphere, where collisional excitations decrease with density, 
while photoionisation remains efficient. We checked this last 
argument by comparing the average ionisation rates from low 
levels (below 2/3 of the ionisation potential) with those from 
higher levels, but the final word will come from the detailed 
analysis of ionisation times, level by level, and using a detailed 
method like the one of Proffitt et al. ( 1999), which will be done 
in future work. 

Clearly, we are confronted with a serious dilemma: which 
method to choose when they yield such different results with- 
out a definitive argument in favour of one or the other? It seems 
reasonable to consider that the MM method gives an upper limit 
of the radiative acceleration; however, when it becomes larger 
than the acceleration given by the light-GLAM method, the cor- 
rect acceleration should lie somewhere in between. To compute 
the diffusion velocities used in the next subsections, we thus de- 
cided to use g™ d values given by the minimum of the respective 
MM acceleration and of the geometrical mean of the accelera- 
tion obtained by both methods: 

gf = Min ^ M ygZ? AMJ ) (18) 

where g r ^ Mi and g r ^ AMi stand for the radiative accelerations 
obtained by the respective methods. 

5.2. Estimates of stratification trends 

The respective rates of abundance increase for each element 
are presented in Figs.|2]and[3] The curves show the logarithm 
of ao(cgs) x 3.1556 • 10 7 (sees) which is the abundance change 
in one year (oo is given by Eq. H2\ ). We have plotted this quan- 
tity only for positive values of ao- Interpreting the curves phys- 
ically, one could say that they give the amount by which the 
abundance of the element considered is increased in one year 
of continuous diffusion. Implicitly we assume that the acceler- 
ations remain constant during this time since we use oo f° r a 
linear extrapolation of the abundance changes by simply mul- 
tiplying with At. This is of course a rough approximation; it 
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is known that the accelerations decrease non-linearly with in- 
creasing abundances when lines are saturated. Furthermore, we 
have not considered any boundary conditions here as has to 
be done when solving the continuity equation. Therefore the 
curves presented here merely display initial trends in abun- 
dance increases. In other words, when the plotted rate is equal 
to unity at a given depth, this means that the abundance of el- 
ement A is doubled at that depth in one year of diffusion. One 
could also consider Figs.[5]and[3]as giving an approximate pic- 
ture of the accumulation of the various elements in the atmo- 
sphere after one year of diffusion. 

Note that if the flux derivative is positive, the trend in abun- 
dance changes can be positive even in places where elements 
are not supported. This corresponds to the accumulation of 
falling particles. Possibly, that kind of accumulation will not 
survive in any fully time-dependent calculation. 

The holes correspond to places where the abundance de- 
creases. Assuming that the abundance stratification process 
tends to a stationary solution (a x — » 0), which still has to be 
demonstrated, a much better prediction of element accumula- 
tions can be obtained through an iterative procedure where the 
accelerations are recomputed after each iteration with updated 
abundance values (see Hui Bon Hoa et al. 2002); this is, how- 
ever, outside the scope of the present work. 

Figure|5] shows the rates of abundance changes for 2 dif- 
ferent values of a horizontal magnetic field (10000G and 
20000 G), as well as the zero field case (angle-independent). 
For some elements, the rates of abundance changes seem to be 
completely insensitive to the magnetic field. This is due to the 
fact that these elements, according to our atomic data, experi- 
ence radiative accelerations mainly through b-b absorptions in 
the neutral state. Therefore the diffusion velocity is dominated 
by the neutral term in Eq.0for which / s / OM , o = 1 holds exactly. 
This is, so far, the only systematic behaviour we have been able 
to clearly identify. For other elements, trends in accumulations 
are highly sensitive to magnetic field strength and orientation 
through fsiowj an d in Eqs. Q and As shown in Fig. [3] for a 
given magnetic field of 20000 G, many elements exhibit signif- 
icantly different accumulations at the magnetic poles and equa- 
tors: the thickness of the respective accumulation zones (or ele- 
ment clouds) often appears to be more uniform and larger near 
the equators than near the poles. Keep in mind that the differ- 
ences between the dotted lines of Fig.|5]and Fig.[3]can only be 
due to Zeeman amplification. We note that ao, which depends 
on the derivative of the diffusion fluxes, is not very sensitive to 
the effect of Zeeman amplification on the accelerations. This 
is because Zeeman amplifications are more efficient in chang- 
ing the maximum overabundance of those elements that can 
be supported by the radiation field than in changing the spa- 
tial derivatives of the accelerations (in Eq. (ll2» . Although we 
expect the final amount of accumulations for each element to 
depend much more significantly on Zeeman amplifications, we 
cannot demonstrate this effect at present since, as mentioned 
above, we extrapolate the ao in a linear way. 



5.3. The silicon case 

At this point it is tempting to examine what happens to silicon 
and to compare these first results to the peculiar behaviour of 
Si discussed in Sect. 14. II According to Figs.|2and|3 a cloud of 
Si starts to form above r as 10~ 2 in the magnetic case; in the 
non-magnetic case, Si starts to accumulate between r « 10~ 2 
and r « 10. This is apparently at variance with Vauclair et al. 
(11979b . who predict that Si does not accumulate in the zero 
field case (or in a vertical magnetic field) but will accumulate 
only when a horizontal magnetic field is present. 

To understand this unexpected result we carried out some 
additional computations assuming 10 times more silicon than 
in the solar case (adopting 7.55 for the logarithmic Si solar 
abundance on a scale where the hydrogen abundance is 12). 
Figure^ displays the relevant data for the analysis of the be- 
haviour of Si. Panel (a) shows the radiative accelerations, ac- 
cording to Eqs. Jl 3I > and (1181 . Our radiative accelerations are 
slightly stronger around r = 1 than those of Vauclair et al. 
(119791 (their Fig. 2c) for 0G, and closer to those of Alecian & 
Vauclair 1981 (their Fig. la). These differences are not surpris- 
ing since our atomic data and our radiative transfer calcula- 
tions are much more accurate. As predicted, due to the pres- 
ence of a horizontal magnetic field, more Si can be supported 
than without a magnetic field for r « 10~ 2 , but the radiative 
accelerations remain smaller than gravity for t < 1 . Silicon en- 
hanced 10 times is not at all supported by the radiation field 
throughout the atmosphere at G; despite somewhat higher ac- 
celeration values above t « 10 _1 , this also remains true for 
the lOkG case. Considering the diffusion flux (according to 
Eq. H0\ ) of silicon in panel (b), one can readily verify that the 
fluxes are positive in places where the radiative accelerations 
exceed gravity (logg = 4.0). For the magnetic case, diffusion 
fluxes go to zero at small optical depths. This can be explained 
by the strong decrease, above t « 1, of the average diffusion 
coefficient (D) in the magnetic case with cos 9 = 0, as shown in 
panel (c) - the relative populations of the Si ions are displayed 
in panel (d). The presence of neutrals mitigates this decrease, 
but not enough to prevent the rather strong decrease in (D) due 
to ions. 

One can conclude that Si is better supported above t » 0.1 
by the radiation field in places where the magnetic field lines 
are horizontal, but less clearly than was predicted by Vauclair 
et al. (I1979> . The total radiative acceleration is less enhanced 
by the magnetic field than expected, and our results reveal that 
Si could also be slightly overabundant in places where the mag- 
netic field is vertical (or zero). On the other hand, we have to 
point out the limits of the "ao x Af" linear analysis carried out 
in this paper, since the complex Si behaviour is not revealed in 
Figs.|2]and|3] 

6. Conclusions 

For this study we have significantly increased the amount of 
physical processes included in the CARAT code. We added 
the radiative accelerations due to bound-free transitions by us- 
ing detailed b-f transition cross sections (whenever available), 
we added the computation of detailed ionisation/recombination 
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rates, and took the momentum redistribution among ions into 
account. This has allowed us to obtain results of higher astro- 
physical relevance than just radiative accelerations. We are now 
able to calculate diffusion velocities (and particle fluxes) in the 
presence of magnetic fields at arbitrary angles. Finally we es- 
timated and discussed abundance stratification trends in mag- 
netic Ap stars. 

After our previous study of the magnetic amplification of 
radiative accelerations, this work constitutes a second step to- 
wards the modelling of abundance stratifications in CP stars. 
Several further steps are still necessary to reach that final goal, 
the most critical of them being an improvement in the theory 
of momentum redistribution among ions. The use of Eq. Jl 8i 
is unsatisfactory, at least for elements for which large discrep- 
ancies are found between the light-GLAM and the MM meth- 
ods (see Sect. 15. U . We thus plan to apply a detailed GLAM 
method following Proffitt et al. ( 1999) but adapted to the mag- 
netic case, which will become possible thanks to a forthcoming 
NLTE version of CARAT. In order to study tepid chemically 
peculiar stars where neutral hydrogen cannot be neglected, we 
plan to implement ambipolar diffusion in CARAT. A good esti- 
mate of final abundance stratifications requires more than sim- 
ple linear extrapolation of diffusion fluxes, since radiative ac- 
celerations depend non-linearly on chemical abundances and 
on the element fluxes at the boundaries (ingoing and outgoing 
particles at the top/bottom of the atmosphere). Assuming that 
a unique stationary solution exists in static atmospheres, an it- 
erative method will converge towards that solution, but in the 
absence of a unique and/or stationary solution, only the simul- 
taneous time-dependent solution of the continuity equation will 
allow us to find the correct solution. 
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Appendix A: The redistribution of momentum, the 
"MM" method 

The collision time for ions can be defined as the time required 
for a deflection by n/2 through collisions with the other par- 
ticles in the medium (mainly protons). This collision time is 
related to the diffusion coefficient (Spitzer 1968 1 by: 

Aifn„ 

U = —pf- Di, (A.l) 
kT 

where D, is the diffusion coefficient of the ion i (presented in 
more detail in Sect.0}, A,- m p the mass (cgs) of the ion, and T, k 
are the temperature and the Boltzmann constant, respectively. 

Let {3 co ij be the collision rate defined as Let f3 ion i and 
/3 rec j be the respective ionisation and recombination rates (see 
AppendixO- The total interaction rate /3, ,j is then defined as 

ficolj "t" Pionj Prec,i' 
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According to the redistribution model proposed by MM, 
and limiting the redistribution of momentum to only the ad- 
jacent ionisation stages (we used the expressions proposed by 
Alecian & Vauclair ll983> . the average redistributed radiative 
acceleration of A° and A +l (for i > 0) may be written as: 



rad Pcolfi 

8o ~ gbb.o + g A + 



N\ Prec,\ Pcolfi 



rad 



Pn 

gA+Ci-i) ■ 

+ 8a^- 



No Plot, l Ptotfi 



Ni-i Pion,i-\ Pcol.i 



Ni Ptot,i-\ Ptot,i 
Ni+l Prec,M PcolJ 



+ gA* 



PcoiJ 
PtotJ 



Ni Ptot,i+\ Ptot,i 

where g A+ > = g b b,i + gbfj (see Sect.|TT 



(A.2) 



(A.3) 



the parameter g h we follow Seaton (2005, private communica- 
tion) and assume go = 0.07, gi = 0.25 (for i > 0). Again, de- 
tailed balancing arguments can be applied to obtain the down- 
ward rate as 



C u = (Ni/NjYdk, 

where the asterisk denotes LTE values. 

The ionisation rate of AppendixlAlis then given by: 



PionJ = JJ- Yj Nik ( Rik + Cik 



ft) 



(C4) 



(C5) 



The recombination rate has been evaluated assuming LTE 
as in Hui-Bon-Hoa et al. J1996i . 



Appendix B: The redistribution of momentum, the 
"light-GLAM" method 

In this method discussed in Sect. 12.21 one defines an energy seg- 
regation threshold such that all the momentum gained through 
a b-b transition of A +l is kept by the ion when the energy of the 
higher level is lower than a given value (here 2/3 of the ioni- 
sation potential). We denote the corresponding acceleration by 
g'hbr The momentum gained through absorption from higher 
energy levels is entirely given to the ion A +(,+1) , and the corre- 
sponding acceleration is denoted by g" hh ,-. The average redis- 
tributed radiative acceleration of A and A +! (for i > 0) using 
the light-GLAM method may then be written as: 

gf^g'bb.^ (B.l) 

g^^^±g" bb ._ 1+g ' bb . +gbfti . (B.2) 



Appendix C: The ionisation and recombination 
rates 

The direct photoionisation rate for ion A + ', in level k, per atom 
is (Mihalas fT978l : 



Rik =4tt I 



hv 



(CI) 



where a v denotes the photoionisation cross section, J v the mean 
intensity of radiation, h the Planck constant, and vo the fre- 
quency of the photoionisation edge. In LTE, J v corresponds to 
the Planck function B v . 

The radiative recombination rate per atom can be derived 
from detailed balancing arguments and we write 

Ru = 4 n \ ± (2hv 3 /c 2 + J v ) e~ hv ' kT dv. (C.2) 

Jv " v 

Approximate collisional ionisation rates based on photoionisa- 
tion cross sections are given by 

C ik = 1.55 \Q n N e T- 1/2 gi a(v ) exp(- M „)/ M „ (C.3) 



where N e is the electron density, or(vo) the threshold photoion- 
isation cross section, and uq = hvo/kT (see Mihalas 1978 1. For 




Fig. 2. Logarithm of the rates of abundance increase (per year) vs. logarithm of the optical depth (at 5000 A) for a horizontal 
magnetic field (90°). All plots have the same scale as the one in the bottom left corner. Solid lines pertain to 20000 G, dashed 
lines to 10000 G, and dotted lines to G. 




Fig. 3. Same as Fig.|2]but for a magnetic field of 20000 G with various inclinations, viz. 90° (solid lines), 60° (dashed lines), 
and ° (dotted lines). 
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Fig. 4. Detailed results for silicon. The dashed lines pertain to OG, the solid lines to 10000G (horizontal). The thin lines are for 
solar silicon abundance, the heavy lines for a 1 dex overabundance (panels (a) and (b)). Panel (a): Radiative acceleration of Si 
(Eq. ( H8» . Panel (b): Diffusion flux of Si (Eq. dlOt ). Panel (c): Average diffusion coefficient {D) = £ NifsiowjD, / N for horizontal 
magnetic fields (cos 6 = 0). Panel (d): Relative populations of the first five Si ions. 



